Accepted for publication in The Astrophysical Journal Letters 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



THE TRANSITING EXOPLANET HOST STAR GJ 436: A TEST OF STELLAR EVOLUTION MODELS IN 
THE LOWER MAIN SEQUENCE, AND REVISED PLANETARY PARAMETERS 
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ABSTRACT 



Cambridge, MA 02138 



Knowledge of the stellar parameters for the parent stars of transiting exoplanets is pre-requisite for 
establishing the planet properties themselves, and often relies on stellar evolution models. GJ 436, 
which is orbited by a transiting Neptune-mass object, presents a difficult case because it is an M 
dwarf. Stellar models in this mass regime are not as reliable as for higher mass stars, and tend to 
underestimate the radius. Here we use constraints from published transit light curve solutions for 
GJ 436 along with other spectroscopic quantities to show how the models can still be used to infer 
the mass and radius accurately, and at the same time allow the radius discrepancy to be estimated. 
Similar systems should be found during the upcoming Kepler mission, and could provide in this way 
valuable constraints to stellar evolution models in the lower main sequence. The stellar mass and 
radius of GJ 436 are A/* = 0.452tE5:oi2 M© and i?* = 0.464t[] ™? R©, and the radius is 10% larger 
than predicted by the standard models, in agreement with previous results from well studied double- 
lined eclipsing binaries. We obtain an improved planet mass and radius of Mp = 23.17 ± 0.79 M© 

and an orbital semimajor axis of 



and Rp = 4.22;°;?^; R, 



a density of pp = 1.69to'i2 g cm 

a = 0.02872 ± 0.00027 AU. 

Subject headings: planetary systems — stars: evolution — stars: fundamental parameters 
individual (GJ 436) — stars: low-mass, brown dwarfs 



stars: 



1. INTRODUCTION 

The applications of stellar evolution theory to astro- 
physics are so widespread, and its validity so often taken 
for granted, that it is easy to forget that it took decades 
to develop, and significant effort to validate by compari- 
son with careful measurements, a process that still con- 
tinues. It is usually only when those theoretical pre- 
dictions fail that the classical discipline of stellar evolu- 
tion "makes the headlines", and even then it draws the 
attention of relatively few. One such instance has oc- 
curred for low mass main-sequence stars. Over the last 
10 years or so it has become clear that our understand- 
ing of the structure and evolution of these objects is still 
incomplete. Discrepancies between theory and observa- 
tion in the r a dii of s tars u nder 1 M©, first mentioned by 
iHoxic (1973[). lLacvl (|1977f) . and others, are now well doc- 
umented for s everal low-mass eclipsing binaries (see, e.g. . 
Popper 19 971: IClausen et al.lll999l: iTorres fc Ribasl l2002t 
Ribas, .2003riL6pez-Morales fc RibasI I2005D . Differences 



in the effective temperatures have been observed as well. 
The direction of these disagreements is such that model 
radii are underestimated by roughly 10%, while effective 
temperatures are overestimated. 

In recent years stellar evolution has had important 
applications in the field of transiting extrasolar plan- 
ets. This is because the planetary parameters of interest 
(mass Mp, radius Rp) depend rather directly on those of 
the star (M^, i?*), and in most cases models provide the 
only means of determining the latter. The s ubject of this 
aper is G J 436, a late-type star found by iButler et al.l 
20041) to be orbited by a Neptune-mass planet with a 
pe riod of 2.644 days. This object was later discovered 
bv lGillon et al.l (|2007^ to undergo transits, enabling its 
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size to be determined (^4 R©). As the only M dwarf 
among the 22 currently known transiting planet host 
stars, GJ 436 (M2.5V) presents a special challenge for 
establishing the stellar parameters, because of the dis- 
agreements noted above. Little mention seems to have 
been made of this, and for the most part past studies 
have relied instead on empirical mass-luminosity [M—L) 
relations to set the mass of GJ 436. Radius estimates 
have often rested on the assumption of numerical equal- 
ity between M^, and i?* for M stars. Despite being the 
closest transiting planet system (only 10 pc away), it is 
rather surprisi ng that the mass of the star is only know n 
to about 10% (|Maness et al.ll2007t [Gillon et al.ll2007bD . 
This is currently limiting the precision of the planetary 
mass, and some of that uncertainty translates also to the 
radius. These two properties are critical for studying the 
structure of the object. 

Given the importance of GJ 436 as the parent star of 
the only Neptune-mass transiting exoplanet found so far, 
and hence the closest analog to our Earth with a mass 
and radius determination, one of the motivations of this 
paper is to improve the precision of the stellar and plan- 
etary parameters by making use of additional observa- 
tional constraints not used before. Specifically, we in- 
corporate the information on the stellar d ensity directly 
availa ble from the transit light curve ([Sozzetti et al.l 
|2007() , which provides a strong handle on the size of the 
star. 

The nature of the discrepancies between evolutionary 
models and observations for low-mass stars has been ex- 
amined recently from both the observational and the- 
oretical points of vie w by iLopez-MoralesI ()2007t ) and 
IChabrier et al.l (|2007t l. Further progress depends on 
gathering more evidence to supplement the few available 
highly accurate mass and radius measurements based on 
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double-lined eclipsing binaries. Thus, a second motiva- 
tion for this work, despite the fact that GJ 436 is not a 
double-lined eclipsing binary, is to present a way of using 
all observational constraints simultaneously to show that 
the star presents the same radius anomaly found for the 
other systems, or more generally, to test the models. Be- 
cause similar constraints may become available in the fu- 
ture for other M dwarfs given the keen interest in finding 
smaller and smaller transiting planets, we anticipate that 
the indirect technique described here may yield valuable 
information on this problem and eventually help improve 
our understanding of low-mass stars. 

2. CONSTRAINTS AND METHODOLOGY 

The basic procedure for establishing the mass and ra- 
dius of a star that is not in a double-lined eclipsing bi- 
nary is to place it on an H-R diagram using observational 
constraints, and compare it with s tellar evolution mod - 
els. We adopt here the models bv iBaraffe et al.l (|1998[ ) 
for a mixing- length parameter of aML — 1-0, which are 
widely used for low-mass stars. The constraints avail- 
able for GJ 4 3 6 are several. The spectroscopic study by 
iManess et al.l (|2007t ) established the effective tempera- 
ture to be Teff = 3350 ± 300 K, and the m e tallici ty was 
estimated photometrically bv lBonfils et al.l (l^OOS") to be 
[Fe/H] — —0.03 ± 0.20. The well-determined Hipparcos 
parallax is tthip — 97.73±2.27mas. ^ yith available visual 
and n ear infrared photometry from iLeggett fc Hawkins! 
()1988l ) and 2MASS, the absolute magnitudes become 
My = 10.610 ± 0.051 and Mk = 6.048 ± 0.052. Ad- 
ditionally, in view of the difRculties in dete rmining effec- 
tive t emperatures for M dwarfs (see, e.g.. iManess et al.1 
l2007t) . we consider as we ll the infrared co lor J — K = 
0.802 ± 0.024. Following [GaxESnterl (pO?)! the 2MASS 
m agnitudes have be en transformed to the CIT s ystem 
of lElias et all ()1982[ ) adopted in the IBaraffe et all (tl998f) 
models, and averaged with those of ILeggett fc Hawkind 
()1988f l. already on that system. The agreement between 
the two sources is excellent. 

Transit light cur ves for GJ 436 have been obtained 
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2007aD in the V band 
20071 ) and lGillon et al l 



from the ground b ylGillon et al.l 

and als o at 8/Ltm bv lDeming et al.l , 

(l2007bD using the Spitzer Space Telescope. Aside from 
minor corrections due to limb-darkening, transit light 
curves can in general be described using three param- 
eters: the radius ratio between the planet and the star 
(Rp/Ri,), the normalized planet-star separation [a/Ri,), 
and the impact parameter {h = acosi/Rj,), where i is the 
inclina tion angle of the orbit. ISeager fc Mallen-Ornelaj 
(|2003f) have shown that a/i?* is directly related to the 
density of the star, and thus c ontai ns valuable informa- 
tion on its size. ISozzetti et al.l ()2007l ) have described how 
a/Ri, can be used together with Teff and stellar evolution 
models to infer M^, and i?*. Briefly, the measured val- 
ues of a/Ri, and Toff are compared with a fine grid of 
model isochrones for a wide range of ages and metallici- 
ties. Theoretical stellar properties are interpolated along 
each isochrone using a small step in mass, and all points 
in the H-R diagram matching the observations within 
their uncertainties are recorded. The best-fitting mass 
and corresponding radius are assigned errors based on 
the full range of model values that are consistent with 
the observations. Other stellar properties can then be 
read off the best fitting model. This is the procedure 
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Fig. 1. — Observational constraints for GJ 436 provided by a/R^, 
Teff, Mfc, and J—K , shown in four different combinations against 
IBaraffe et al.l II1998I) isochrones for solar composition. Ages be- 
tween 1 and 10 Gyr are represented, although the curves are hardly 
distinguishable because stars in this mass range evolve very slowly. 
The raggedness of the isochrones in the right panels is due to the 
limited precision with which magnitudes are tabulated in the pub- 
lished models. The agreement between the latter and the observa- 
tions appears good in all panels, but the implied masses and radii 
show large systematic differences indicative of internal inconsisten- 
cies in the models. 

we apply below. For GJ 436 we have restricted the 
comparis on to solar metallic ity, given the [Fe/H] esti- 
mate by iBonfils et all (j2005D . The values of the light 
curve parameters we adopt are weighted averages of the 
results from the ground-based and Spitzer photometry: 
a/R^ = 13.34 ± 0.58, Rp/R^ = 0.0834 ± 0.0007, and 
6 = 0.848 ±0.010. 

3. MASS AND RADIUS DETERMINATIONS 

In previous studies the mass and radius of GJ 436 have 
been derived in three ways: either 1) the mass has been 
obtained fro m the near-infr ared ( JHK) mass-luminosity 
relations of Delfossc et al. (2000) (giving = 0.44 ± 
0.04 M0: IManess et al...2007.) and the radius has been as- 
sumed to be numerically equal to the mass (jGillon et al.l 
l2007al) . or 2) the mass has been held fixed at the above 
value and the radius cons trained directly from the light 
curve ("Gi llon et al.l l2007a.b). or 3) both Af* and R^ have 
been solved for simultaneously subject to the constraint 
that they be numerically equal and making use of the 
implicit sensitivity R cx M^^^ between mass and radius 
in the light curve fi tting procedure ([Deming et al. 2003; 
Gillo n et aLli2007bD . In the second case the radius ob- 
tained is near 0.46 R© for M^, = 0.44 Mq, and in the 
third case = i?^ « 0.47 or 0.48 in solar units. Thus, 
differences of order 0.02 to 0.04 remain between these 
determinations, depending on the procedure used. 

The methodology in the present paper is completely 
different, and can yield improved precision and also give 
a better understanding of possible systematics. We ini- 
tially applied the a/Rj, and Teff constraints as described 
in fj2|using the lBaraffe et al.l (1998) models, and obtained 
a mass value near = 0.50 M©, which is consider- 
ably larger than previous estimates. The predicted abso- 
lute visual magnitude (My = 9.86) is also much brighter 
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than that computed directly from the Hipparcos paral- 
lax. This inconsistency strongly suggests a problem with 
the models, which is not entirely unexpected for a star 
of this type. As an alternative to a/i?*, we then exper- 
imented using the absolute K magnitude as a proxy for 
luminosity, as well as replacing Toff with the color index 
J—K} Figured] displays the four constraints in different 
combinations against solar-metallicity isochrones from 1 
to 10 Gyr. The seemingly good fit in all planes belies the 
serious discrepancies present in other derived quantities 
that are not shown explicitly. Those results are listed in 
Table [H Masses inferred from a/i?* are systematically 
larger than those from Mk^ and so are the radii. The 
masses from Mk come close to the estimates from the 
empirical M—L relations, but the corresponding radii are 
considerably smaller than expected. This would seem to 
go in the direction of the results from eclipsing binaries 
(see ©. 

On the other hand, there is good evidence from vari- 
ous sources that the bolometric luminosities from these 
models are not seriously in error ( Delfosse et al. 20001 : 
[Torres fc Ribai 120021: iRibasI 120061 : iTorres et al.l I2OO6I ). 
This suggests that adjustments to the model radii and 
temperatures might resolve the discrepancies in Table [U 
and allow us to obtain a meaningful result for GJ 436. 
We explored this by introducing a correction factor /? to 
the radii, and at the same time applying a factor /3~^/^ 
to the temperatures in order to preserve the bolometric 
luminosity. We repeated the comparison between the ad- 
justed models and each of the four sets of constraints for 
a range of /3 factors centered on the value indicated by 
the eclipsing binary studies. Figure [5] shows the result 
for several of the key stellar properties. The lines corre- 
sponding to the four sets of constraints seem to converge 
for a value of (3 near 1.1 (representing a 10% correc- 
tion to the model radii), which happens to be the typical 
factor found by the eclipsing binary studies mentioned 
earlier. For this value of fi the models yield essentially 
the same mass, radius, and luminosity for GJ 436, in- 
dependently of which set of observational constraints is 
used, as one would expect from a realistic model. Thus, a 
self-consistent solution is achieved. To arrive at the best 
possible values of Mi, and i?* we next applied all the con- 
straints simultaneously, and varied /3 as before, seeking 
the best agreement with the measurements. The result is 
illustrated in the bottom right panel of Figure [U where 
the quality of the match as represented by is shown 
as a function of the correction factor over the restricted 
range in which the models agree with all four observables 
within their errors. The best match is again near /3 — 1.1. 
The resulting mass and radius are M^, = 0.452l;g Qi[2 ^& 
and i?^ = QAMtolfi R©. These and other inferred stel- 
lar properties for GJ 436 are listed in the top section 
of Table [2] We emphasize that these quantities are the 
result of the simultaneous application of the four con- 
straints, and the agreement with some of the values in 
Table [1] is accidental. 

^ We have refrained from directly applying any constraint based 
on the V magnitude because of suspected deficiencies in the models 
for optical passbands, related t o missing molecular opacity sources 
short ward of 1 ^m (see, e.g., IBaraffe et al.l [19981 : IDelfosse et al.l 
l2000h . Experiments using the I—K color as an alternative indicate 
that the I band is also affected at some level. 
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Fig. 2. — Stellar properties for GJ 4 36 derived from the com- 
parison with stellar evolution models by IBaraffe et all II1998I) . as a 
function of the adjustment factor applied to the model radii (/3). 
The observational constraints given by a/i?*, Tgff , Mjf, and J—K 
are applied and shown in pairwise combinations, as labeled. Kinks 
in the curves are a reflection of the discreteness of some of the 
quantities tabulated in the models. The simultaneous application 
of all four constraints results in the curve shown at the bottom 
right, indicating a best fit for a /3 value near 1.1 (see text). This 
value is represented with the vertical dotted line running through 
this and the other panels. 

4. THE RADIUS DISAGREEMENT WITH THE MODELS 

Our mass and radius estimates are compared with mea- 
surements for late-type double-lined eclipsing binaries in 
Figure [21 Only systems with the most accurate determi- 
nations are shown (relative error s below 3%), which are 
taken from the summary by iRiba s (2006) . For compari - 
son we include two isochrones from .Baraffe et al.l ([1998f ) 
corresponding to ages of 300 Myr (as estimated for two 
of these binaries) and 3 Gyr (more representative of the 
field) . All eclipsing binary systems are seen to have larger 
sizes than predicted for their mass. We note also that 
GJ 436 lies in the gap between masses of 0.43 Mq (for 
CU Cnc A) and 0.60 M© (GU Boo B), and thus provides 
valuable additional information on the radius discrepan- 
ci es for low mass stars. 

[L6pez-Morales[ ()2007l) has investigated how these dis- 
crepancies (Ai?*/i?,t) depend on metallicity and the 
strength of the chromospheric activity, quantified in 
terms of the X-ray luminosity (specifically, Lx/^boi), for 
the rather limited sample available so far. Both of these 
factors have been suggested to play a role. The above 
study examined single M dwarfs as well as M dwarfs in 
binar y systems. GJ 436 is a rather inacti ve star for its 
type (lEndl et al.l[2003[ : [Butler et al.[[2004 ). but was de- 
tected nonetheless by ROSAT as an X-ray source because 
of its proximi t y. Th e X-ray luminosity was reported by 
[Htinsch et all |l999( ) to be Lx = 0.7 x 10^^ erg s'^. 
When combined with the bolometric luminosity in Ta- 
ble d we obtain LxlL\,oi = 7.() x 10" ^. The metaUicity 
was estimated by [Bonfils et al.l (|2005[ ) to be near solar: 
[Fc/H] = -0.03 ± 0.20. Considering GJ 436 as a single 
star, the value Ai?*/i?* ~ 10% we find is consiste nt with 
the overall conclusions of [Lopez-MoralesI ()2007( ) in the 
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Fig. 3. — Mass-radius relation for all double-lined eclipsin g bina- 
ries w ith relative mass and radius errors under 3% (data from lRibasI 
120061) . GJ 436 is shown at the values determined from our model- 
ing (open circle), and is seen to display the same radius discrepancy 
as the other s ystem s. Two solar-metallicity model isochrones by 
IBaraffe et all lll998l ') are shown for reference, corresponding to ages 
as labeled. 

sense that it is similar to the offsets for other systems 
regardless of ix/^boi, and at the same time it seems to 
follow the trend with [Fe/H] exhibited by other single M 
dwarfs. 

5. PLANET PARAMETERS AND FINAL REMARKS 

With the host star properties known, the planet 
parameters we infer are given in the bottom section 
of Table [2l The required orbital period and vcloc- 
ity se mi-amplitude K^, are adopted from Mancss ct al. 



( 2007 ). along with the eccentricity from Demor v et al.l 
2007| ). who obtain a similar value for Ki,. The im- 



proved precision of these derived parameters is a re- 
flection of the better stellar parameters. The slightly 
larger planet radius than in previous studies confirms 
with even greater statistical significance the conclusions 
of earlier authors re garding the presence of a hydro- 
gen/h elium envelo pe ( Gillon et al.l 1200763 : iDeming et al.l 
I2OO7I: iCTillon et al.ll2007bD and agrees very well with the 
models bv iFortnev et all ((20071) for a 10% fraction of 
those elements. 

In this paper we have shown that GJ 436 provides a 
valuable test of stellar evolution theory near the bottom 



of the main sequence, made possible by the fact that it 
has a transiting planet. Traditional studies in the area of 
low-mass stars have made the comparison with models 
by measuring the mass and radius directly for double- 
lined eclipsing systems containing M dwarfs. In a few 
other cases angular diameters have been measured in- 
terferometrically for single stars, and the ma ss has been 
inferred from empirical M — L relations fe.g.. iLane et al.l 
l2001t IS^gransan et al.ll2003l ). More recently, a variety of 
constraints and assumptions have been used to infer the 
mass and radius of the late-type secondaries in F-l-M 
syste ms observed as part of transiting p l anet surveys 
(e^ lBouchv et all 120051 : iPont et ani2005t iBeatty et all 
|2007| ). Though perhaps not as compelling as having ac- 
tual model-independent mass and radius measurements, 
the approach in the present work is able to make use of 
available information for GJ 436 and compare the models 
directly with the observational constraints without re- 
quiring a direct measurement of the mass and radius. 
The discrepancy in i?* is derived by parameterizing it in 
terms of a single adjustment factor to the model radii 
assuming the luminosity from theory is accurate, as 
other observations seem to indicate. 

NASA's upcoming Kepler mission, currently slated to 
launch in early 2009, will emphasize the search for tran- 
siting Earth-size planets. These should be easier to de- 
tect around late-type stars. Therefore, we anticipate that 
many systems similar to GJ 436 could be found and 
become a significant source of information on radii for 
low-mass stars, since they will have all the observational 
constraints needed (including trigonometric parallaxes) 
to test models of stellar evolution in the way we have 
done here. 
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TABLE 1 

Stellar parameters for GJ 436 ba sed on different sets of constraints, using the 
MODELS by IBaraffe et al] (jl998l '). 



Observational constraints 



Parameter Teff and a/Ri, J — K and a/Ri, Tf.g and Mk J — K and Mk 



M* (Mq) 


407+0.003 


n c:it;+0.042 
U.010_Q 


448+0-018 


452+0-014 


R* (R©) 






421+0 009 
'^•^^1-0.009 


421+0-009 


log 9* (cgs) 


. 702+0 020 


4.780t°«g 


4 841+0-020 
^•"^1-0. 009 


4 843+0-018 


a/R^ 






14 63+0-34 


14 65+0-31 

l'±.00_g 21 


Tcff (K) 




3684l«^ 




3585^23 


J — K (mag)^ 


• 0.801«:«g? 




n sin+0-002 




Mv (mag) 




q 71+0.27 
'^-0.36 


in 256+0-0™ 

iU.^OD_Q 1594 


10 244+0-082 


Mk {raa-g)^ 


r- yco+0.090 


a.DD_Q 22 






(L©) 


n n'?4«+o.oooo 


n n-^s-^+o.oogg 

U.UOOO_g QQgg 


0260+0-0015 


0260+0-001-1 



^ Magnitudes are in the CIT photometric system of lElias eFall 11983) . 



TABLE 2 

Stellar and planetary parameters for the 
gj 436 system. 



Parameter Value 

Stellar parameters 

(Mq) 0.45210-014 

(R©) 0.46410-0?? 

i*(L0) 0.0260j:noi7 

logff* (cgs) 4.843l°;0}f 

Age(Gyr)'^ 6l* 

Planetary parameters 

Mp (Me) 23.17 ± 0.79 

(R©) 4.2210:0? 

Pp(gcm-3) l-69lo.}2 

logffp (cgs) 3.107 ± 0.040 

a (AU) 0.02872 ± 0.00027 

^ Due to the unevolved nature of GJ 436 the age is 
essentially unconstrained by the observations. We list 
this value only for completeness. 
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